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["T i i Abstract: The present status of the theoretical description of Initial (ISI) 

and Final State Interactions (FSI) in nucleon induced exclusive meson pro- 
duction processes close to threshold is reviewed. Shortcomings are ad- 
■ dressed and new perspectives given. 

> 

O 

oo : 1 Motivation 

O 

It is quite surprising, that the discussion of the theoretical description of 
ISI and FSI in nuclear and particle collision/production processes is gaining 
again increasing attention, even though people thought, this issue has been 
well understood due to an established formalism going back to outstanding 
scientists such as e.g. K. Brueckner, G. Chew, E. Hart, K.M. Watson, A.B. 
Migdal, E. Fermi, G. Gamov and A. Sommerfeld entering nearly all stan- 
dard text books of Scattering Theory, e.g. P3E]- The existing formalism for 
the treatment of ISI and FSI is based on the idea, that T-matrix elements 
of many processes (in the Distorted Wave Born Approximation (DWBA)) 
are well described as overlap integrals of initial and final state wavefunctions 
multiplied by some interaction potential. This idea works quite successfully 
for (mainly nonrelativistic) textbook examples, in which excitation ener- 
gies are small compared to the masses of the particles involved, i.e. for 
the description of electromagnetic transitions in atoms (excitation energies 
~ 1 eV) or /^-decays/transitions of atomic nuclei (excitation energies ~ 
10 2 . . .10 3 keV). In these processes the initial/final states hardly go offshell 
due to the transition interaction, or at least for a very short time. In re- 
action processes like NN — ► NNX (with X being a mesonic system) close 

1 e-mail: kleefeld@cfif.ist.utl.pt, URL: |http://cfif.i st. utl.pt/~kleefeld/ 

2 Invited contribution to the symposium on Threshold meson production in pp and pd 
interaction, extended COSY-11 collaboration meeting, 20.-24.6.2001, Cracow, Poland. 



1 



to the particle production threshold the situation is quite different. The 
excitation energy of the NN system in such a process is of the order of the 
mass of X. In the case of reproduction this yields already about half of 
the nucleon mass, for heavier mesons the excitation energy is in the range 
of the nucleon mass. At high excitation energies is a difficult task for the 
outgoing NNX system to rearrange the 3-momenta and bring the 3 outgo- 
ing particles to their mass-shells. In the case of an incoming or outgoing pp 
system the infinite range of the Coulomb force even prevents the initial or 
final state particles respectively to go to their mass shell within a finite time. 
The production of "light" 7r-mesons seems to be quite exceptional. Yet even 
there one can see already various new aspects being important for the de- 
scription of ISI and FSI of processes under consideration. The issue of this 
presentation is not the discussion of various experimental signals for mainly 
FSI in NN — * NNX shown e.g. in |S] (and references therein). It is 
a introductory review of present theoretical formalisms and new theoretical 
developments treating ISI/FSI to enable researchers involved to find to a 
(still outstanding) quantitative covariant theoretical formulation to describe 
ISI/FSI and a better quantitative analysis of experimental data. 



2 Invariants and notation in NN — ► NNX 

In the following I will introduce quantities used in the theoretical descrip- 
tion of the reaction NN — > NNX. For convenience I use the metric 
(+,—,—,—), denote the outgoing meson X by <fi, the incoming/outgoing 
particles by 1, 2, 1', 2' and 3', i.e. JV(1)JV(2) -> N(1')N(2') 0(3'), and the 
respective 4- momenta by pi, p%, Pv, P2' and p-y = k^. 5 independent 
Lorentz invariants s, si, S2, ti, ti and 1 dependent invariant S3 are in- 
troduced by s := (pi +P2) 2 , si := (py + P2') 2 , s 2 := (p 2 ' + P3>) 2 , h ■= 
(pi ~ Pv) 2 , t 2 := ( P 2 ~ Pv) 2 and s 3 := (p 3 , + p v f = m\, + m| + m\ - 
s\ — S2- For a discussion of the 3-body phasespace I refer the reader e.g. 
to 00 IE]- Using the kinematical triangle function \(x,y,z) := x 2 + y 2 + 
z 2 — 2 (xy + yz + zx) and defining l^lcm := <J A(s, m^, si)/(4 s) there are 
various ways to introduce quantities, which carry the same energy informa- 
tion as the cm-energy y/s, e.g. Q cm := y/s — my — my — m^, \pi\ cm = 
\p2\cm ■= y / A(s,mf,m |)/(4s), k := v /\(s,m 'l,m 2 .)/(2(m 1 + m 2 )) and 
V* ■= \h\™ X /™<t> = ^X(s,m 2 ,sr n )/(ism 2 ) with s™" = {m v + my) 2 . 
Here kq is the Lab-wavenumber of the relative motion of the incoming N N- 
system. The excess-energy Q cm and the dimensionless quantity which 
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is the maximum possible cm-momentum of the produced meson in units of 
the meson mass ms, vanish exactly at threshold. The relative motion of the 
outgoing 2-particle subsystems 1'2', 1'3' and 2'3' can be described e.g. in 
terms of relative Lab-wavenumbers K23, K31 and k := Kyi'- 

y / A(si,mf,,m|) ^A(s 2 ,m|,,mp JA(s 3 ,mf,,m|) 

'~ 2(m 1 >+m 2 >) ' K23 ' 2(m2' + m^,) ' K31 ' 2(m 1 '+m ( j ) ) 

(1) 

3 The Watson-Migdal Approach 

Historical remark 

The present theoretical formalism of treating ISI and FSI goes mainly back 
to ideas formulated by K. Brueckner, G. Chew and E. Hart (1951) [S], 
K.M. Watson (1951,1952) [13 E], A.B. Migdal (1955) QU and E. Fermi 
(1955) ^21) an d has been transported through the literature by now as 
Watson-Migdal Approximation (WMA) or Distorted Wave Born Approxi- 
mation (DWBA). The idea behind this formalism is, that the T-matrix Tfi 
of a particle production/anihilation/scattering process can be "factorised" 
into bound-state or (elastic) scattering wavefunctions describing the parti- 
cles in the initial and/or final state, and a short-ranged interaction part T^ ' , 
describing the transition from the initial to the final state. The T-matrix in 
the Watson-Migdal Approach can therefore be factorised in dimensionless 
enhancement factors Tpsi and Tisi describing the influence of FSI and ISI 
respectively (being 1 for no FSI and ISI respectively) and a short-ranged T- 
matrix describing the production/scattering process without ISI/FSI: 

Tfi a TfsiT^Tjsi. (2) 

In Watson's first extensive approach JJ| considering a process 12 — ► 1'2' the 
energy dependence of the T-matrix and therefore of the (differential) cross 
section da/dn stems mainly from the 2-particle final state wavefunction ipfLj 
which is asymptotically related to the wavenumber k, phaseshifts 5l and the 
scattering length ol (L = orbital angular momentum of the l'2' system) 3 : 

e t5L sm.5 L a L k l 

l fi <XTPf L ~ kL+1 - K 2L+l( cot(5i _ i ) " l + ia L K 2 W { ) 



3 K. Brueckner et al. developed the idea, that Tfsi is given by the wavefunctions at 
the interaction point r — 0, i.e. by ipf(0). This idea has been later used by R. van Royen 
and V.F. Weisskopf |71) in the initial state to estimate decay rates T by T oc |-!/jj(0)| 2 . 
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Here I used the Effective Range Expansion (ERE) k 2L+1 cot 5l = — a^ 1 + 
0(n 2 ) ^1] . As a consequence Watson obtained for the square of the T-matrix 
and the differential cross section the following wavenumber dependences: 

2 \a L \ 2 K 2L da sin 2 5 L 2 2 \a L \ 2 k 2L+2 

l J M ^ 1 + | aL |2 K 4L+2 ^ « K 2L «" l J /il « 1 + | aL |2 K 4L+2 " 

(4) 

Watson's arguments have been generalised to processes with 3 particles in 
the final state. Let's e.g. consider the process NN — ► NN<f> close to thresh- 
old. Assuming for simplicity that FSI are mainly described by the outgoing 
destorted NN wavefunction ijjf,L, while representing the produced meson 
4> by a spherical wave, which is for the orbital angular momentum I of the 
meson relative to the NN system just proportional to a spherical Bessel 

function je(q' r) 9 oc q a (with q' := |fc^| cm ), one could — according to 
Watson — suggest Tfi oc q' e if)f^L- Including the 3-body phasespace one 
obtains within a nonrelativistic framework for the total cross section: 



a oc 



r m <t>'n<t> rm^rjA, I 1 2 2 L+l 



Using the threshold identity « ^ (m 2 rj 2 - q' 2 ) ( m *^™ N + 0(q' T 

a Taylor expansion in m yields for the cross section close to threshold: 

a^NN^ [[NN] L <p] e ) oc V 2 / +L) + 4 . (6) 

Derivation of the Onshell-Watson-Migdal Approximation 

Now I am going to derive a more sophisticated Watson-Migdal formalism, 
which yields the traditional Watson-Migdal Approach in a certain limit 4 . 
The derivation — already fully performed by Watson — has recently been 
rederived and critically reconsidered by the author in |151 1161 117] . while 
simultaneously promoted in ^H]. Let's consider a system of colliding and 
produced particles being described by a Hamilton operator H = K + V = 
K + W + U , with K being the kinetic energy operator, W being the particle 
number nonconserving short range interaction potential and U being the 
particle number conserving interaction potential of long range. Furthermore 

4 The formalism to be presented is essentially equivalent to the so called two potential 
scattering formalism which can be found in various textbooks of quantum scattering theory 
(see e.g. Throughout such a formalism the potential consists of a sum of a short- and 
long-ranged potential. 
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I define the long range Hamilton operator h = K + U. The respective 
eigenstates to the operators H, h and K fulfil the following Schrodinger 
equations: 

(E - H)\iP± > = , (E- h)\xt > = , (E n - K)\<p n > = (7) 
yielding the Lippmann-Schwinger equations (G^(E,K) := (E — K zizie)^ 1 ): 

\ti>= \<p*>+ E _B ±ie v\<p a >= \xi> + E _ l H±l w\xi> 

\xt>= \<Pa> + „ ... U\<p a >= \ip a >+G±(E,K)T±(E,K)\<p a > 
Hi — n rt %E 

The S-matrix elements of the considered process are determined by Sp a = 

< > yielding the T-matrix Tp a = < ^~\(U + W)\^ > . For E = 
Ei = Ef the T-matrix can be reformulated in terms of Watson's Theorem 
(see e.g. p. 448, 461 in 0), i.e.: 

Tp a = < Xp\U\^ a > + < x*\W\il>+ > = < <Pf)\U\x+ > + < il>a\W\xZ > 

(8) 

As the particle number conserving potential U is positioned in matrix el- 
ements of states with different particle number, the respective matrix el- 
ements have to vanish. As a result I obtain Tp a = < xjl^lV'a > = 

< V^l^lxi >• I n the following step we perform the first approxima- 
tion widely known as Watson-Migdal Approximation (WMA) 5 . The exact 
asymptotic states \tp^ > are approximated by the [7-destorted asymptotic 
states \xt >; i- e - IV'a >— \Xa >■ As a consequence I obtain an approxi- 
mated Watson's theorem Tp a ~< \p~ \W\x% >• This approximation is valid 
to first order in W and to all orders in U. Defining the free propagator 
G^(E, K) = (E-K± ie)- 1 and inserting complete sets of free asymptotic 
states \<j) n > the T-matrix in the WMA can be separated into a short- and 
long-ranged part via: 



2> a ~ <<p \(l + Y l T£(E,K)\<p n >G+(E,E n )<<p n \) W 

n 

'l + Yl > G o( E > E ™) < ^l T e!(£' K )) \<Pa > • (9) 



5 This approximation yields a "factorisable" T-matrix for a field-theory, which is by 
construction known to be not factorisable. To the knowledge of the author there is no 
theoretical work to quantify the validity of this approximation. Even the consequences 
of the approximation to the unitarity of the approximated S-matrix for systems of many 
coupled channels have to be investigated. 
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If the incoming/outgoing particles are only shortly /weekly offshell, one may 
approximate the free propagator Gq(E,K) in the following way 6 : 



G ° ±(jB ' K) = E-K±ie = P E~K T ™ S ( E - K ) ra Ti*S(E-K). 

(10) 

The remaining ^-distribution enforces the incoming or outgoing particles 
respectively to go to their mass shell and to scatter elastically. This is 
why I say, the remaining principle value part of the propagator contains 
all the "offshell" information of the scattering/production process, which is 
in problems with many coupled channels needed to restore the unitarity of 
the S-matrix 7 . Let's see, how the OWMA shows up in a non-relativistic 
treatment of a simple 12 — * 1'2' scattering process, in which both the initial 
and final state are treated onshell. In order to proceed I identify the free 
asymptotic states \(p n > with momentum eigenstates \Hi > and |«/ > in the 
initial and final state respectively, while in the intermediate state I insert 
complete sets of momentum eigenstates |<& > and \q*f > with J d 3 qi \% >< %\ 
= f d?qf\qf >< <jf\ = 1 resembling the symbolic completeness relation 
\<Pn >< (fn\ = 1- Using the nonrelativistic dispersion relation E = 
n't /(2m) = nf /(2m) the T-matrix Q in the OWMA is given explicitly by: 



, _ _ . onshell 
T(Kf,Ki) ~ 

r I ~* |2 

K (<X f \-i«j d\<K f \T+(E,K)\q f > 6(E-fL) <q f \)w 

/I -> |2 
d 3 q t \qi > 5{E - M-) < %\T+(E,K)\Ki > ) . (11) 



6 This approximation I will call here Onshell-Watson-Migdal Approximation (OWMA). 

7 K. Nakayama is surely right, if he states, that offshell effects have not to be observable 
within a field-theoretical formalism. He claims, that results of the principle value inte- 
grals are cutoff dependent and therefore dependent on the renormalisation scheme. To the 
opinion of the author this is due to the approximations made within the Watson-Migdal 
Approach. A careful parametrisation of the (unknown) principle value contribution of the 
propagator restoring field-theoretical constraints like unitarity of the S-matrix etc. will 
again remove the renormalisation-scheme dependence of the Watson-Migdal Approach. 
The discussion about observability of "offshell" effects shows, how important and instruc- 
tive a careful quantitative formulation of the treatment of IS1 and FSI is. 
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Now one can perform a partial wave expansion of the elastic T-matrices 
in the initial and final state. The result is given in the footnote 8 . It is 
obviously of the form (J2J), i.e. T/j ~ Tpsi Tjsi- It is straight forward 
to generalise the result to the "offshell case" , which contains also information 
about the principle values in the propagators G^(E, K). Parametrising the 
principle value contribution by some unknown function P{n) and denoting 
:= < ^j, Lf\W\Ki, Li > the T-matrix in the WMA can be expressed as 
follows 9 : 

TiR^Lf^U) = Tpsi(Lf) T(°) T ISI {Li) = 

l + ie i5 *t iK,) sm5[ f (K f ) (1 - P f (K f )j) 
\ + ie iS UM sin^^Oa- 
l + i(e 2a 4W-l)(l-P /(K/ )))T(») 
l + ^(e 2i ^ (Kl) - 1) (1 - Piim))) (13) 



The axially symmetric partial wave expansions of the elastic T-matrices are given by: 
-(2^) 2 m <Zf\l$(E,K)\<tf> = 

= 7" E (2L, + l)e V"'' sin 5/ .(«/) *(*/.?/)) 

-(2 7r) 2 m <q i \T+ l (E,K)\FU> = 
1 °° 

= - Y, (2Li + l)e iS ^ Ki} tmi6i t (Ki)P t , i (a» 
Denoting T(k/; Ki) ~ YIl l- T(J$f, Lf', «i, Li) I finally obtain in the O WMA: 



T(Kf,Lf,Ki,Li) — 



= (l + ie^f^ smS[ f {K f )) < K f ,L f \W\^U > (l + ie is ^ Ki) sin 
cot 5/ (k/) cot 



cotSl f (K f )-i J ' Jr ' ' cot «i 4 («,)-< 

^!2!< %w >w|w. (12 , 

Here /£. and («:/) denote the well-known Jost functions (see e.g. of the elastic 
scattering problems in the initial and final state respectively. 

9 The offshell parametrisation P(k) is related to the offshell quantity V(k) introduced 
in EH1 by P(n) = iV(K)/(a L k 2L+1 ) with /t 2L+1 cat 6 l (k) = - a" 1 + 0(k 2 ). 
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or 

T{R f ,L f -K U U) = Tpsi(Lf) T(°) TjsdLi) = 

cot6l f ( Kf ) - iP f ( Kf ) t(o) cot^(«Q - iP,( K ,) 
coiS[(K f )-i cot 5 l.(Ki) -i 

Refl f ( Kf ) + iP f ( Kf )lm fj f (K f ) Refl^+iPij^lm fjjnj) 

//,(*/) 

(14) 

It is interesting to see 2 ways how to "switch off" ISI or FSI: one could 
perform the limit 5l — > 0, or one could set the P(/c) — ► 1. The OWMA is 
obtained by P{k) — * 0. In the limit P(k) — * — 1 the enhancement factors 
just reduce to a multiplicative complex phase 10 , i.e. Tisi — ► exp (2 idl (Ki)) 
or Tps/ — > exp (SiS^ («/))• It has to be explored in future, in how far the 
restoration of the unitarity of the S-matrix will set some restrictions on the 
allowed values of the offshell parametrisation and its limits, and between 
the offshell parametrisations in the initial and final state 11 . Tjsi and Tpsi 
derived above are also valid for complex phaseshifts, i.e. even in the case of 
inelasticities yielding \rji,\ < 1 for rjz '■= exp(2i5i) = \t)l\ exp(2 i Re Sl). 

Enhancement factors, 2-potential formalism and unitarity of 
the S-matrix 

To get a better understanding of the properties of enhancement factors I 
want to "define" them in a very schematic way within the standard 2- 
potential formalism of Quantum Scattering Theory (see e.g. 0). For a 
1-channel scattering problem with a potential V, which is a sum of a short- 

10 The careful reader may be reminded in the so called Fermi-Watson Theorem (see 
e.g. based on the idea of the unitarity of the S-matrix, which claims Tfi = 

exp(iRe<5/) \T/i \ exp(iRe<5i) and Im<5/ =Im<5;. 

11 Among the first few models studying the principle value contribution or offshell ef- 
fects in ISI and FSI there are the following three to be mentioned: Ch. Hanhart and K. 
Nakayama 18 investigate the "half-shell" function related to P(k) within a separable po- 
tential model. V. Baru et al. |19II2U| find essential differences between separable potential 
models and realistic NN potentials. They calculate P(k) in TViV — » NNX for a simple 
rescattering graph. It would be very interesting to see, what will change in their results, 
if the produced meson is coupling to an excited resonance of the nucleon. Finally I want 
to refer to the work of B.O. Kerbikov |21| . He compares for the first time in a systematic 
way the effect of "half-shell" and the Jost functions to enhancement factors. 
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and long-ranged potential, i.e. V = Vs + Vl, the S-matrix S is given in terms 
of the corresponding phaseshifts by S = exp (2 % 5) = exp (2 i (5s+5l)), while 
the T-matrix can be decomposed in the standard manner: 



sin 5 



iSr 



sin 5l 



+ e 



2i5r 



sin 5s 



T L + e 



us. 



--Tr 



(15) 

Let's now consider a many-channel scattering problem. Here the unitary 
S-matrix is determined in terms of "phaseshift matrices" A, As and by 
S = (S" _1 ) + = exp (2 i A) = exp(2i(A 5 + A L )) such, that A = A 5 + A L is 
Hermitian. The question is now, whether it is possible to define an invertible 
matrix Tpsj = Tf SI such, that S = exp (2 i A) = Tpsi exp (2 i As) Tjsi = 
(5 _1 ) + . If As is Hermitian (yielding of course also the Hermiticity of 



leads to Tp SI 



-fsi 



and 



Al), the unitarity of the S-matrix S = (S~ 
Tjgj = Tfg T . Remembering, that due to the Baker-Campell-Hausdorff for- 
mula we can shift the argument x of a function / according to f(x + a) = 
exp (a gj) f(x) exp (-a ^) by a differential operator J| obeying , x] = 1, 
we can ask the question, whether there exist a "functional derivative ma- 



trix" 
S = 



-g^j such, that 
exp f 2 i (As + Al 



8 
6A S 



Ac 



1, yielding a factorisation of the S-matrix 



)\ and T-matrix T = (S — l)/(2in), i.e.: 



5 
T 



exp |^A L 
exp ( A L 



SA S 
5 



Tfsi 



exp ( 2 % As ) exp 



T 5 exp - Al 



5 



5A f 



Tisi 



SA S 



(16) 



It is important to realise that the matrix has to act like a functional 

T fsi en " 



derivative 12 . On the other hand the unitarity condition T F l 



forces 



5 

5A S 



A, 



FSI 



0. We can learn from the "naive" considerations pre- 



sented above 3 points 



13. 



1 One easily can proof, that there exists for a finite dimensional matrix As no matrix 
M such, that the commutator [M , As] yields the unit matrix. 

13 It is left to future work to investigate the constraints of analyticity properties of 5- 
particle amplitudes like the amplitude of the process NN — * NNX to enhancement 
factors. For the discussion of 2- and 3-body unitarity in elastic T- and S-matrices I refer 
the reader e.g. to |2^1l2l] . 
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• It is a very difficult — most probably only approximately soluble — task 
to systematically construct enhancement factors, which fulfil simultaneously 



the requirements 



8 

8A S 



Ac 



1 and 



8 

<5A S 



A, 



0. The first commutator 



can be fulfilled only on a functional or operator level. 

• Tpsi and Tjsi are related (Tpsi = Tfgj) and constrained by unitarity (see 
also [TZHHI). 

• The construction of enhancement factors is strongly based on A<j, i.e. 
enhancement factors do depend on short-range interactions 141 . 



4 Energy expansion and behaviour of the total 
cross section at close to threshold 

General idea 

The total cross section of NN — > NN<f> is given by (defining P := p\ + P2, 
P' := pv + P2> + and the flux factor T{s) := 2 \J A(s; m'f, m 2 ,)): 

S f d 3 pi> d 3 p 2 > d 3 k 

a{s) ~ J(s) J (2vr)3 2^(|^|) (2vr) 3 2^(|p r |)( 27 r)3 2^(|^|) 

• (2vr) 4 5\P'-P) \ff^ ■ (17) 

The symmetry factor 5 is ^ for two identical nucleons in the final state, 
otherwise 1. |7/j| 2 denotes e.g. spin averaging, i.e. \ Ylsi s 2 s, s, l^7«| 2 - 
As for this reaction o~(s) at threshold vanishes exactly 15 it is suitable to 
perform an energy expansion of the a(s) in a sense of a Taylor expansion in 
terms of the small energy dependent parameters 7u or Q cm . The argument 
of the phasespace integration \Tfi\ 2 = \Tpsi Tisi\ 2 can be expressed 

14 With respect to the last point, I have to admit, that V. Baru et al. |19U20| are right, 
if they state |20) : "The absolute magnitude of FSI effects depends on the momentum 
transfer (or on the mass of the produced meson) and hence is not universal." They 
deduce from a simple rescattering graph in iViV — > NNX, that the enhancement factor 
only for heavy produced mesons seem to be universally factorisable. Yet if they state, 
that ". . . FSI effects cannot be factorised from the production amplitude . they ignore 
the possibility, that for every short-ranged sub-amplitude in the process NN —* NNX 
the enhancement factor may be constructed individually in a matrix form discussed. Of 
course it seems more attractive to treat the process NN — > NNX in a Faddeev formalism 
rather than factorising short- and long-ranged interactions. Yet for processes involving 
more particles (especially with spin) it seems at present beyond the numerical capabilites 
to avoid a Watson-Migdal formalism. 

15 This is not true for nucleon nucleus collisions, which allow also subthreshold produc- 
tion cross section due to e.g. the Fermi motion of the nucleons in the nucleus. 
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as a function of the 5 independent Lorentz- invariants s, s%, s 2 , t\, t 2 . As 
discussed above, close to threshold Tpsi will mainly depend on n\ 2 , K23, K31, 
i.e. Tps/ ~ T FS i(s,si,s 2 ). T IS i depends mainly on k , i.e. T/ 5 / ~ T IS i(s). 
If the propagators in the short-ranged T-matrix T^°'(s, s\, s 2 , t\, t 2 ) are not 
too close to their pole and particles are mainly produced through nucleons 
or S-wave resonances, it is a good approximation to treat close to threshold 
the T-matrix as a constant 16 , i.e. T<® ~ T^(s thr , s{ hr , s% hr , t{ hr , tf T ) with 

s thr = (my + m 2 < + m^) 2 , s\ hr = (my + m 2 >) 2 , = (my + m^) 2 , 
frf = m\ — my (m 2 i + m^), t^ r = m\ — m^ (my + m 2 /) . (18) 

To improve this approximation one can perform a Taylor expansion of 
T(°)(s, si, s 2 , ti, t 2 ) in the variables 5s := s — s thr , 5si := s\ — sf 17 " ', 5s 2 := 
s 2 — s|' ir , St\ := t\ — tl hr , 8t 2 := t 2 — t 2 ir , or simply a Taylor expansion in 
the variable 5s := s — s thr , i.e. 17 : 

TW(s, Si , 53,^2) = 

= t ^ m 2^ tUt2> = f»" T^,s 2 , tu t2 ,, 19 , 

n=0 n=0 

Instead of performing a Taylor expansion in 5S or Q cm it is more suitable 
to perform a Taylor expansion in the dimensionless variable t?^, which is 
related to s and Q cm by (2/i := (mi' + m 2 >)/m ( j } ): 




Q 



^ +1 + 2^ = ^ 2 + l + ^ + (2^. (20) 

Pulling all factors depending only on s in front of the phasespace integration 
and selecting one specific spin-isospin channel one obtains for cr(s): 

1 9 00 

i 2u N (\py\) 2 u N (\p 2 ,\) 2^(1^1) 
■ \T F si(s,s 1 ,s 2 )T^(s 1 ,s 2 ,t 1 ,t 2 )\ 2 . (21) 



16 This approximation yielding very compact analytical results has been shown to work 
quite well for the reactions pN — » pNX with X — rj, tt° (for details see e.g. 1 151 172| 1. A 
corresponding approximation has been later applied to similar reactions by I25| . 

17 This kind of a energy expansion ot is (in terms of the variable r]$) for the first 
time taken into account by |26j . The problem is, that the Taylor coefficients are derived 
from onshell meson nucleon scattering data and not from any theoretical meson nucleon 
potential. 
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The energy behaviour of a(s) is determined on one hand by the s-dependent 
function in front of the phasespace integral consisting of the flux factor 
the initial state enhancement factor Tjsi and the factor (5s) n , on the other 
hand by the energy dependence of the phasespace integral itself, which is 
affected strongly by Tpgl- I n the following I perform the threshold energy 
expansion of the various factors in terms of the varibable rj^ separately. 
The flux factor behaves very close to threshold like a constant, while the 
dominant term in the power expansion of 5s is the term with (<5s)°, as all 
terms with n > will be of higher order in 77^. 

Treatment of ISI 

The energy expansion of Tjsi(s) = 1 + 2 ( e J<5 ( K o) sm <5(k ))(1 — -P(«o)) is dis- 
cussed in |16j . I will give here merely the general ideas. Even if many people 
think, the offshell effects in the initial state are very small, i.e. P(kq) — 0, 
the relation Tpsi = Tf$j shows, that this is most probably not true. There- 
fore I perform the Taylor expansion P(kq) = P^ + P^ rj\ + 0(7/4). As 
phaseshifts 18 are given in terms of the kinetic Lab-energy T^ab = (s — (m\ + 
7B2) 2 )/(2m2) (particle 2 at rest!), one performs the following Taylor expan- 
sion of the data around threshold: 

S(T Lab ) = 5(Tl^)+d'(TZ b )(T Lab -T^l) + 0((T Lah -T%Z) 2 ) 

Ynp" u / 1 \ ^ 

= S (Tffi+6'(l$S,)-±Z (1 + — ) ijJ+Ofo}) 

= 5^+5^4+0(4). (22) 

C^Lab = (( m l' + m 2' + m <t>) 2 ~ ( m i + m2) 2 )/(2 m2)). Combining all we get: 

T IS i{s) = 1 +i(e i5W sin<y(°))(l- P (0) ) + 
+ i (e 2is(0) 6W(l- P(°))- pV)e i5i0) sinjW) 4+ O(^) . (23) 

The NN phaseshifts at the threshold of pp — > ppX are badly known for 
mesons heavier than the pion. The accuracy of the determination of high 
energy phaseshifts limits a quantitative theoretical calculation of cr(s). ISI 
- close to threshold — don't seem to affect the energy dependence of the 
cross section. Their effect is just a multiplication of the cross section by 
some number 19 (see also e.g. |28p. 

18 I treat here phaseshifts as complex numbers containing also eventual inelasticities! 

19 Assuming P(k ) ~ and taking S (0) w -60° and S (2) w from the VPI-phaseshift 
analysis |27| one obtaines for pp — > ppr\ close to threshold \Tisi(s )\ ~ 1/2, i.e. a 
reduction of a(s) by a factor of about 1/4. 
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Treatment of FSI 



From considerations presented above I conclude, that in leading order the 
energy-, i.e. 77^-dependence of a(s) close to threshold is determined by 
the energy-, i.e. the ^-dependence of the phasespace integral over a FSI- 
dependent integrand. I.e. the problem is now to determine the threshold 
energy- or 77^-dependence of the following FSI-modified phasespace integral: 

R fsi (s) . = f d3 Pi' d3 P2' d*k S 4 (pl _ p) . 
3 W J 2^(1^1)2^(1^1)2^(1^1) V ' 

• \T FS i{s,s 1 ,s 2 )Ti°\s 1 ,s 2 ,t l ,t 2 )\ 2 . (24) 

The phasespace integration can be rewritten into a restricted integration 
over four Lorentz-invariants sx, s 2 , ti, t 2 , or into a restricted integration 
over three of these invariants (e.g. si, s 2 , tx) and one unrestricted angle, 
while the fourth of these invariants (e.g. t 2 ) is expressed in terms of the three 
and the angle (see e.g. [§]). Remembering, that \Tn\s\, s 2 , tx, t 2 )\ close to 
threshold is nearly constant, I want to investigate now only the threshold 
energy dependence of the phasespace integral with a integrand / being a 
function of sx, s 2 and s only 20 , i.e. (S := s/rrii, Sx := si/m 2 , S 2 := s 2 /m^): 



R pir s) „ f d*Pv d*pv d*k 6 i {p ,_ p) f ( S ,Sx,S 2 ) . 
3 {) 7 2^(1^1)2^(1^1)2^(1^1) 1 >JK,l,2) 

(25) 

The case considered most in NN — ► NNX close to threshold is the sce- 
nario, when Tpsi is described as a sum of T-matrices describing the FSI 
of the 2-particle subsystems 1/2', 2'3' and 1'3' in the final state only, i.e. 
Tfsi(1'2'3') ~ Tfsi(1'2') + T FSI (2'3') + T FSI (3'1') - 2. This scenario in 
the OWMA is equivalent of using only the lowest order terms in the Fad- 
deev expansion of the 3-particle elastic T-matrix in the final state. In an 
incoherent language 21 with the correct limit for no FSI this would yield 
I Wl'2'3')! 2 ^ I Wl'2')| 2 + I W2'3')| 2 + I W3'l')| 2 - 2. Hence in 
this case we have / (S, Sx, S 2 ) ~ fi 2 (n 12 ) + fiz(^2z) + /3i(«3l), i-e.: 

RFSifs) ~ f d3pv d3p2 ' d3k 8HP' P) 
3 W J 2^,(1^1)2^(1^1)2^(1^1) V ' 

• (/l 2 (« 12 ) + /23(^2 3 ) + /3l(«3l)) • (26) 

20 The investigation of the general case is devoted to future research. 
21 The coherent case for my 7^ m.2' has still to be investigated. 
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If the functions /#(/%) can be Taylor expanded, i.e. fij(Kij) = Yl c ij K ij> 

a 

then the threshold behaviour of R^ SI (s) is given by |17| : 
7T 2 m ? +2 



R[ SI (s) ~ Y — — r-y — - 
3 v ; ^ 2a+1 ^ 



/ T \ a+4 / t \ a+4\ 

» , ,(«) fm 4 ,J+\ ( a ) (m <j> J-\ \ f„<*H,nf„<*-w 



(27) 

using the definitions A := y 7 ]. — (mi/ — m,2') 2 (mi/ + m 2 ')~ 2 and J Q := 

J^du ((1 + 2/i)A 2 (l and J± := ± ^2 ± 2 VT^A* + (2/j)A 2 . 

Considering these results it is easy to see 22 , that the free 3-body phasespace 
(i.e. f(S, Si, S2) = 1) is proportional to 77^. Finally I want to discuss the 
more general case (|25|l. yet for equal masses my = my = m^, i.e. A = 1. 
In this case R^ SI (s) can be rewritten as |15j : 

= J—n 2 ml V lf l du ® / dv \ [f{S,Si,A + Bv) 

2VS Jo J r] 2 u + 1 J z 

+ f(S,Si, A-Bv)](28) 

with VS = y/r% + l + ^ + (2/i) 2 , fiTi = 5+1- V45 ^u + l, A := 

i [5-5i+2/j 2 + 1] and 5 := ^5r/ 2 u (5i - (2/i) 2 )/5i. Remembering that 
5, Si, ^4 and B are functions of 77^ and u, i.e. S 1 = S(r]V), Si = Si(u,rjV), 
A = A(u, 77?) and i? = B(u,r] 2 ), it is easy to see, that the integrand can 
be Taylor expanded in the variable 77^. Afterwards one has to calculate 
analytically/numerically the u-/-u-integrals in the Taylor coefficients. Using 
Mathematica^ it is no problem handle Taylor expansion and integrations. 

The Coulomb problem and how to surround it 

Upto now we learned, that we easily can deduce the threshold 77^-dependence 
of the total cross section by performing a Taylor expansion of \Tpsi T^\ 2 in 
the relative wavenumbers, say e.g. k. This is always possible for well behaved 

22 Using fgdu \Ju(l — u) = 7r/8 one easily can deduce the expression for the free 
phasespace, i.e. R 3 (s) = fi m 2 A r\% + 0(r,%) . 
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long-ranged potentials, which yield e.g. for the S-wave the standard ERE of 
the following form (see e.g. p. 466 in [23], [3U] or p. 12ff in |31|): 

K cot 8 (k) = - - + - r k 2 - P r 3 k 4 + Q r 5 K e + ... (29) 

If we keep only the leading two terms in the ERE, i.e. only k cot 8o(k) = 
— a -1 + r k 2 /2 we call this ERE "shape independent" . Specific properties of 
the long- and / or short-ranged potential under consideration can absorbed in 
higher order "shape dependent" terms of the expansion. In 1959 M. Cini, S. 
Fubini, A. Stanghellini (CFS) |32| included the singularity of an One-Pion 
Exchange (OPE) potential at k 2 = -ml /A to the ERE by (see also 

Kcot 5 (k) = + \r k 2 - - P K — 7T . (30) 

a 2 1 + q k 1 

The absorption of more pathologic singularities like the essential singular- 
ity of the one-photon exchange (at k 2 = 0) in the Coulomb interaction is 
performed within a formalism using a "Modified Effective Range Function" 
(MERF) Kff( K 2 ) := A l (k 2 ) + B l (k 2 ) k 2L+1 cot 8 l (k) [331 dp. The cor- 
responding "Modified Effective Range Expansion (MERE)" is given by: 

B l (k 2 ) k 2L+1 cot 5 l {k) + A l (k 2 ) = -- + \rK 2 + ... (31) 

It is easy to see, how the MERE enters Tpsi described e.g. in (|14jl. i.e.: 

K 2 L+1 (cot 5 L (k) -iP(n)) 



Tfsi 



K 2 L+ 1 (cot 5 L (k) -i) 

- ~a + § « 2 - • • • - A M - i k 2L+1 B l (k) P(k) 



(32) 



So, if \Tpsi\ 2 can be Taylor expanded, we immediately can derive the thresh- 
old energy dependence of the total cross section. To handle e.g. the pp-FSl in 
PP ~^ ppft , we have to use a Coulomb Modified Effective Range Expansion 
(CMERE) suggested by G. Breit et al. (1936) and extensively investi- 
gated e.g. by |13 E3 EE1 EH 24 . For S-wave scattering the CMERE yields 
Bq(k) = Cq(?]) := 2 7r ?7/(exp(2 7r rf) — 1) ("Coulomb penetration factor", 
"Gamov factor") and Aq(k) = 2Kr/H(r]) with the "Sommerfeld parameter" 
r\ := a/vpp = a \J (k 2 + (m^m^' / (m\> + my)) 2 )/ k 2 and a ~ 1/137 25 . In 

23 As shown in 221 the functions Al(k 2 ) and Bl(k 2 ) are related to the Jost functions 
and Jost solutions of the potential problem under consideration. 

24 For the theory of pp-scattering see e.g. also gU EE), P- 49ff in EU,EZII331- For 
the quantum-mechanical Coulomb problem and the ERE see e.g. [481 1491 I5T3] , p. 133ff in 
|2,E1E1], p. 323 in ED, p. 285 in gU, p. 54ff in [3T], P- 161ff in p. 16 in [Tul].) 

25 Note, that H(r]) = -7-ln»7 + E^Li V 2 /{s (s 2 + V 2 )) and 2kt] ~am v = (28.82 frn) -1 ! 
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the 2-potential formalism the total phaseshift <5o( K ) can be decomposed into 
a sum of the Coulomb phaseshift o" (k) (exp(2z<7o) = r(l + irf)/T(l — irj)) 
and the phaseshift S c (k) of the short-ranged potential in the presence of the 
Coulomb potential, i.e. 5o(k) = ao(n) + 8q(k). In the CFS-CMERE derived 
by D.Y. Wong et al. (1964) jSHj, which is based on a 2-potential formalism, 
one obtains (see e.g. also p. 468 in |22|, p. 1869 in [H], or [111 EUJ El 01113]): 



C 2 ( V ) KC ot5 c ( K ) + 2ktjH( V ) = -- c + - \r c n 2 - f ^ 2 - (33) 

a c 2 1 + q c k z 



The construction of an enhancement factor described e.g. in (|14|) is straight- 
forward: 



Tfsi 



l + -.e 



2 i 8q(k) 



!-)--[ g2jo- (K) g 2 «<5 c (k) 



1 + I f e 2i<r («) — ^ COt ^) + l l 
2 \ k (cot (5 c (k) — i) 



I ) ( ! - P(k) 
1 1 ( I - P(«) 

0(1 /' 




2 i cto(k) 



+ ... -2^7? #(??) + inC 2 (ri) 
+ ... -2kt]H(t]) - !kC 2 (?j) 



1 - P(k 

(34) 



Due to essential singularities appearing in the differentiation of (7q(k) and 
Cq(t]) in the limit k — > it is for a Coulomb potential not possible to 
perform a Taylor expansion of \Tpsi\ 2 in the variable k, in order to determine 
the energy behaviour of the total cross section close to threshold. Yet the 
problem can be surrounded in the following way: first one regularises the 
Coulomb potential according to (fJ,,£ > 0)(g := am p = lim 2kt]): 

K — >0 



U(r) = g 



- /ir 



r + e 



dt 6(t - fj) e 



p.i. 



' g dt [6(t-fi)-e6(t-n)] e^- t)e 



- 1 r 



(35) 
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As the potential fulfils | J* °° drrU(r)\ < oo and is of Yukawa type, one can 
determine by the method described in the corresponding Jost solution 26 
fo(n,r) and Jost function /o(k) := fo(K, 0). As the potential is analytic at 
r = 0, it is then possible — assuming k to be real positive — to apply the 
method in section IV of jSU to determine from fo(K,r) and /o(k) a MERF 



Kq^k 2 ) for the regularised potential by (/q(k, 0) := lim <9/o(k, r)/dr): 

r— »0 



+ ■ 



1 



/o'M) 

/<>(«, 0) ' |/ («)| 



K (cot£jf(«)- 



1 1 2 

j T K + . 

a 2 



(37) 



The regularised enhancement factor is now obviously obtained from (|34[) by 
the replacements e 2 " 70 ^ -» / *(k*)// (k), C 2 (r/) -> 1/|/ (k)| 2 , 2 kt? F(r?) -» 
/o'MV/oM) - WI/o(^)| 2 , 5 C («) - <^M. The result is: 



Tfsi 



/o(«) 



+ ... 



/o'M) 
/o(«,0) 



+ 



2 i K 



+ 



/o(«.Q) 

AM) 



1 1-P(«) 



(38) 



Remember, that / *(«)//o(«) = l/o(^)| 2 /(/o(M) 2 - Now |7> 5/ | 2 can be Tay- 
lor expanded and the threshold energy dependendence of the total cross 
section be determined. The Coulomb FSI is obtained by [i — > +0 and 
e — ► +0. Eq. (|38|) can be also used, to construct enhancement factors from 
Jost functions obtained for singular potentials by the WKB approach [57] . 



Effective Range Expansion and astrophysics 

One valuable source jSHl EH] for sophisticated Coulomb destorted wavefunc- 
tions at low energies and CMEREs to be used in pp scattering is astrophysics, 

26 The Jost solution fo(K,r) is obtained to be (Ei(— z) := — J°° exp(— t)/t for z > 0): 

fo(n,r) = 

exp(- r n (r + e)) 



e IKr {l + ^ (g exp(^e)) n / dr n - 

.Inn K' 



+ f— 1 K + T n 



n 



K 2 + (— i K + T v ) 2 



} 



{i + E ( a ?." e) )i Ei( ~ |r + e) "»> - «"""" + " Ei <-(-- + «)(»« - 2i «»] 

n=l \ 1 / 



n 



7 — z \ +i (- - arctan — ) 



} 



(36) 
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which is due to a interesting correspondence between NN — * NNX close to 
threshold and theory of radiative capture processes in astrophysics 27 . 



5 |7fst| onshell beyond a threshold expansion 
The FSI model of Faldt and Wilkin 

A particularly interesting existing nonrelativistic theoretical model for the 
description of FSI is the model developed by Faldt and Wilkin |6U1 1611 
1621 I63j . It does not only predict the leading terms in the ^-Taylor ex- 
pansion of the total cross section close to threshold, it gives expressions 
describing the threshold energy behaviour of the total cross section to all 
orders in t]a. The idea of the model is the following: suppose we know 
the threshold energy dependence of the differential cross section of a pro- 
cess N1N2 — ► (Ni/N 2 ')bs </>> while (NyN^BS denotes the bound state 
(BS) of two components Ny and A^' with binding energy eb- Then we 
can derive — by analytical continuation of the final state wavefunction 
- the energy dependence of the total cross section of a process N\N 2 — * 
(NyN 2 ')ss 4> with a scattering state (Ny N 2 ')ss of same quantum numbers 
as the BS. According to Faldt and Wilkin the BS-wavefunction and the 
corresponding SS- wavefunction are related through an analytical continua- 
tion by lim {v/2 a (a 2 + k 2 ) i/jss(k, r)} = — ipBs(r) with a = J\xyi £B and 

^12 = m y 1 + m 2' ■ As in a Watson-Migdal picture the T-matrix factorises 
into wavefunctions, there is now the hope, that also the T-matrix elements 



Compare e.g. the reaction pp — » dir + close to threshold and the capture process 
pA — » (A + 1)7. In the first process the final state d-ir + is nearly at rest, while in 
the second process the initial state pA is nearly at rest. The capture process behaves 
like time reversed compared to a threshold production process. Both processes have 
enhancement factors. In NN — > NNX we have Tfsi, in the capture process we have the 
astrophysical S-factor S(E). The T-matrix elements look very similar: for iViV — > NNn 
we have Tfi ~ J dr u; Ji(|fcir| r)in with Uf at E ~ 0, while in the capture process we have 
Tf i ~ f drufr x m with m at E ~ 0. The only difference is, that the excitation energy of 
the meson production process is high, while in the capture process all particles stay nearly 
onshell. 



18 



of the two processes are related by 28 [SH] (</ := l^lcm): 

\TfiiNxN2 -» (N v N 2t ) ss <l>) | 2 \TfiiNtN2 -» (N^N^nst] 



(2tt) 3 2 mi' (2tt) 3 2 m 2 ' 2 a (a 2 + k 2 ) (2tt) 3 2m BS 



(39) 



da(N!N 2 -> (N v Nv) B s4>)/dn oc ^ IT^N^ -> (N v N 2 >)bs <t>)\ 2 « g ,n 
at threshold yields iT^AqiVa -> (N v N v )ss<l>) ? « q ,nr - 1 /(2a{oc l + k 2 )). 
After performing the phasespace integration one obtains for the total cross 
section with the S'S'-system in the final state (/i -1 := (m^ +m2') _1 + ^I 1 ): 



a(JViJV a -> (NyN^sst) oc t# P {n) (v<pCi>) with Q, := • (40) 

Some of the functions P^^rj^Q are listed without and/or with Coulomb- 
corrections in [On]- As mentioned in |63| I17j the cross section <r(pp — ► ppn ) 
close to threshold can be well described by a(pp — > pp7r°) oc rj^P^ 1 '^ £„■) = 
r)\ C|/(4 (1 + yl + ^ 2 C 2 ) 2 ) indicating a quasi-bound state in the correspond- 
ing {pn}/ = i final state. For the discussion of the relation between the 
{pn}/ = o system and the deuteron in the final state see e.g. [HH1 E2J 165] . 
The A^AX-final state is addressed e.g. in |66j . 



Onshell Faddeev models 

For completeness I want to mention here some works going in the direc- 
tion of a full Faddeev calculation for NN -> NNX. Moalem et al. (1995) 
|H7) developed a semiquantitative formalism for describing Tfsi (and 
expanding the elastic 3-body T-matrix in the final state according to a 
Faddeev expansion upto second order. First solutions of simplistic nonrela- 
tivistic Faddeev models using separable interaction kernels can be found for 
the ^-system in |68| 169] and for the reaction pp — * pprj in [201 ■ The main 
result of the works is an investigation of the effective range parameters in 
the r]N system and the rjNN coupling constant. The models have to be 
improved to be quantitative and predictive 29 . 

28 The approximation made in the approach of G. Faldt and C. Wilkin is, that one has to 
match the T-matrix element of a 2 — + 2 process (described by 2 independent the Lorentz- 
invariants s, t) to a T-matrix element of a 2 — » 3 process (described by 5 independent the 
Lorentz-invariants s,Si,S2,ti,t2)- This is only possible by imposing some approximative 
(onshell) constraints to 8,si,sn,ti and/or ti. 

29 It e.g. should be remarked, that the complex r]N scattering length can't be investi- 
gated without considering also the effective range. For the description of a meson nu- 
cleon interaction potential parameters, which reproduce phaseshifts, have to be choosen 
consistently. Finally separable interaction kernels may introduce spurious states (also 
quasi-bound states) not existing in the original theory yet seen in the phaseshifts. 
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6 Relativistic OWMA — the ABSSM 



As discussed in Section the Watson- Migdal Approach has been based on 
the idea, that the energy dependence of the T-matrix close to the reaction 
threshold is mainly driven by ( "factorised" ) destorted wavefunctions of inter- 
acting subsystems in the initial or final state, which are usually determined in 
nonrelativistic Schrddinger frameworks. With the exception of lightfront ap- 
proaches the consistent factorisation of relativistic wavefunctions — so called 
n-particle Bethe-Salpeter Amplitudes — of bound and scattering states from 
the T-matrix has upto now been an open problem 30 . In order to remove this 
lack of formalism and taking into account the fact, that production processes 
of heavy mesonic systems close to threshold involve high energy and momen- 
tum transfers, the author recently developed the relativistic so called Asymp- 
totic Bethe Salpeter State Method (ABSSM) [13 E3 ESI , which has 

been ap- 
plied [151 174j as a first test in a simplistic way to the reaction pp — > p A K + 
in the quark-gluon picture 31 . In a second step the ABSSM is now applied 
to reactions like NN — ► dX. How is the relativistic deuteron wavefunction 
factorised within the ABSSM? A composite state vector \P,B> being nor- 
malised according to <C P',B'\P,B >= (2vr) 3 2oj b (\P |) 5 3 {P'-P) 5 BIB can 
be expanded into free Fock states. The overlap of the free n-particle Fock 
states with \P,B > is then reexpressed in terms of respective n-particle 
Bethe-Salpeter Amplitudes. For details I refer to |73| . The most dominant 
term of the composite deuteron state vector is the two nucleon sector, i.e. 32 : 

\d + (P,M)> ~ 



- EE 

Sl,t\ S2,t2 



[ f P* [ d 3 Xl ^3 e ~i (ft • X! + p 2 ■ X 2 ) 

J (27r) 3 2m JV (2vr) 3 2m JV J 

b + b + (p 2 ,S 2 ,t 2 )\0 > (pi,sx,tx) u (2) (p 2 ,s 2 ,t 2 ) 

«0|T[V (2) (x 2 ) V (1) (xi)]|d+(P,M) > n n . (41) 



30 In decay processes there is the general belief, that due to R. van Royen and 
V.F. Weisskopf |71) the decay rate is — to a good approximation — proportional to 
the square of the wavefunction of the decaying system at its origin. Yet this assumes the 
factorisability of this square of the wavefunction. 

31 It should be mentioned that publication |74| contains apart from a quark-gluon de- 
scription of pp — *■ p AK + the first theoretical investigation of this reaction in the nucleon- 
meson picture taking into account beside the pseudoscalar meson-exchanges (ir, K) the 
upto now by different authors neglected, yet — as shown in [741 — important vector- 
meson exchanges (p,K*). 

32 Here Sj,t, denote spin and isospin. Spinors/operators fulfil u(p 7 s,t) u (p, s' ,t') — 
2mS ss , 5„,, {b(p,s,t),b+(p',s',t')} = (2itf2Lj(\p\)S 3 (p-p>)5 as ,5 tt ,, ... 
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In this form the composite state vector can be conveniently used for the 
calculation of T-matrix elements with the help of Wick's Theorem in the 
standard manner. In phasespace integrations the composite deuteron ap- 
pears only as one particle with mass m^. It is beyond the scope of this 
presentation to discuss all the promising results and relativistic corrections 
predicted by this new method 33 . 

7 What can go wrong? 

The unitarity violation problem 

Naive use of enhancement factors may violate unitarity, especially in many 
channel problems! Tpsi and Tisi are related by unitarity! 

The fitting problem 

Accurate experimental measurements of the total cross section of NN — > 
NNX close to threshold are supposed to prove or disprove theoretical de- 
scriptions of the meson nucleon dynamics. As is seems, that the energy 
dependence of total cross sections close to threshold is mainly described by 
FSI parametrised e.g. by a sophisticated Effective Range Theory taking into 
account eventual quasi-bound states, the experimental measurement of to- 
tal cross section constrain essentially one number, i.e. |Tv°) Tjsi\- Under the 
assumption, that we can estimate \Tisi\ quantitatively, we can constrain 
the modulous of the short-ranged T-matrix |T^|. It is therefore curious 
to observe, that the majority of present theoretical models for the descrip- 
tion of NN — ► NNX introduce or keep one free parameter in their models, 
which allows them to bring their calculated cross sections into the experi- 
mental data. The most distributed method (see e.g. [2U I2H1 EE] , the seventh 
equation in [21] and the second equation in 18 ) of introducing a "fitting 
parameter" is the following: As people know from the WMA, that Tpsi is 
proportional to the dimensional quantity ipf(r = 0) ~ l/(«(cot S — i)), they 
introduce a dimensional constant or function N ("numerator") to obtain 
a dimensionless enhancement factor Tpsi by Tpsi = N/(K(cotS — i)) = 
N/(— a~ l + r k 2 /2 + . . . — i k). This numerator N is either fitted to the 
data or chosen without much reasoning 34 to be — a -1 . It is also possible 

33 A further decisive quantitative test of the ABSSM will be its application to the positro- 
nium decay compared to standard Bethe-Salpeter descriptions or the new factorisation 
method based on a dispersion approach presented in |75|. 

34 As it has been shown in Section^ the numerator is given by N = re cot 8 — i reP(re) = 
— 1/q _1 + r re 2 /2 + . . . — i reP(re) with -P(re) staying undetermined. Even in the OWMA, 
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to transfer the quantitative uncertainty in the description of FSI/ISI to a 
"fitting parameter" in the short ranged T-matrix T^°\ While in the model 
of V. Bernard et al. |25| |T^°^| 2 is just chosen to fit to the data, the "fit- 
ting parameter" of the Jiilich-model for the threshold production of pions 
in nucleon nucleon collisions (see |77| and references therein) is the heavy 
meson nucleon (anti-) nucleon coupling constants (uNN, aNN) in the so 
called "pair diagrams", which is tuned to meet the pp — ► ppn total cross 
section. 

The Gell-Mann-Watson-Rosenfeld-Koltun-Reitan et al. problem 

The interpretation of the total cross section of pp — > ppn close to threshold 
has a peculiar history. The the theoretical situation is summarised by H. 
Machner and J. Haidenbauer (1999) [ZHj as follows: "From phase space 
alone the cross section should follow a . . . rf dependence. However, this is 
in complete disagreement with data. Taking into account 'minimal effects' 
from the final NN interaction i.e. using the first term in the ERE 35 . . ., 
leads to an rf dependence Then the authors give a widely and by now 

used expansion of a{pp — ► ppir ) close to threshold, which they call "barrier 
penetration model", i.e. 

a(pp -» ppir ) = Pss r)l + P Ps rj* + (3 Pp rfi . (42) 

This threshold expansion goes back to a phenomenological discussion of M. 
Gell-Mann, K.M. Watson (1954) jZ| and A.H. Rosenfeld (1954) [HO]. Early 
datapoints for rj^ > 0.5 listed in [BJ seemed to support the ^-behaviour of 
a(pp — > ppir ) very close to threshold, while the 2 datapoints at % = 0.660 
and 1.11 [SOI indicated a(pp — ► ppir ) cx rj^ (see also ^Sl) showing the 
onset of Pp-production. This is why theory thankfully acknowledged the 
theoretical model for a(pp — > ppn ) presented by D.S. Koltun and A. Reitan 
(KR1966) [H2]) which resulted in a(pp — > ppir ) cx r]% close to threshold 36 . 
Accurate data for 0.13 < r] w < 0.56 published by H.O. Meyer et al. (1990) 
|84l I85j and A. Bondar et al. (1995) jHS] could be — within the errorbars — 
understood as a final confirmation of the ^-dependence of a(pp — ► ppir ) 

i.e. with P(k) — 0, there is no reason to neglect in N the term r n 2 /2, while keeping it 
finite in the denominator. Such a procedure strongly affects the energy dependence of the 
cross section. 

35 At this point the authors refer to the expression |Tfs/| 2 = ^/(k 2 + [— l/o + r k 2 /2] 2 ). 
36 It should be mentioned, that 1969 there followed a complementary theoretical ap- 
proach by M.E. Schillaci et al. |S2]- 
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close to threshold 37 . Yet latest data (see jSTJ |HH1 and references therein) 
show clearly, that the energy dependence of cr(pp — ► ppir ) very close to 
threshold is much steeper than rj 2 . It seems, that phasespace behaviour, i.e. 
a 77^-dependence, is restored for % < 0.2, while the "r/^-plateau" of a{pp — > 
PP 7T °)/ r lw m the range 0.2 < rj^ < 0.5 is a structure due to a quasi-bound 
state in the pp-system well described within the formalism of G. Faldt and C. 
Wilkin [HOI URESES] (see also Section^}. It is no surprise, that — according 
to the introduction of [SHj — the most theoretical (standard) models of pp — > 
ppTT are more or less sophisticated versions of the work of KR1966, who 
differ mainly in the sophistication of , yet use same FSI and phasespace 
arguments in order to reproduce the ^-dependence of <r(pp — > ppn ) close 
to threshold. The problem is now, that the calculations of KR1966 and 
therefore also the — often quoted — original considerations of M. Gell- 
Mann, K.M. Watson (1954) |7S] and A.H. Rosenfeld (1954) [SU] contain an 
explicit mathematical mistake. To make the mistake more transparent, I 
translate the notation of KR1966 in the language used in this presentation, 
i.e. I perform the replacements M — > m p , [i — ► m n , p' — > k, E p — > T' := 
Tjf\ n = y/ai — 2m p and Ef — > Q cm - Using this notation equations (18) and 
(19) in KR1966 read: 

3 

da _ 2tt 1 v^|^| 2 dp{Q cm ) dp{Q cm ) _ (2m 7r m p )2 / — _ 

dT'~ v 4 ^ bl dT< ' dT> ~ (2^) 4 VWcm > 

(43) 

Making according to KR1966 the approximation |Tf>| 2 ~ |T(0)| 2 /(1 + a 2 k 2 ) 
and using k 2 = m p T' + 0(T' 2 ) we obtain 



2 4 



4u Jo 1 + cr m p i + G>(i z ) dZ 

(44) 

This should be compared to equation (21) in KR1966. Koltun and Reitan 
write: "For Q cm 3> rn~ l a~ 2 we may approximate Eq. (18) and write": 

27T 1 V^lTVnM^ /" Qcm AT' 1 d P(Qcm) 2 , > 

i.e. they obtain a total cross section with a ry 2 .- and not a 77^ -dependence. We 
can see, how the error occurs: although the lower integration limit is 0, they 



37 If one has a more careful look to the data, one can see already in these data very close 
to threshold a deviation of the energy dependence of a(pp — » ppiv ) from the "expected" 
7^-behaviour. 
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allow the replacement 1 + a 2 m p T' — ► a 2 m p T', which is equivalent neglect- 
ing the —a -1 term in the ERE 38 . As easy the error is detectable in KR1966, 
as hidden it is among various features of T^ ' in the publications and com- 
puter codes of present theoretical calculations. One either can recover it by 
trying to reproduce their results or from the statements they make in their 
publications. Many models mentioned in E. Hernandez et al. j90l I91j . espe- 
cially the models of C.J. Horowitz et al. [021, J-A. Niskanen IHSUHl, G.A. 
Miller et al. 7(i and the Jiilich-model (see [77] and refererences therein), 
are set up nearly in the same way as KR1966. As all these works seem 
to reproduce the energy dependence of the direct pion production close to 
threshold predicted in Koltun and Reitan [§2], it seems, that they all use 
the same approximation as Koltun and Reitan [S2] 39 . Ch. Hanhart states in 
Tab. 1.2 of his doctoral thesis [HE]) that an oc rf and on page 13 below eq. 
(1.7), that ^NN^NNni 7 !) x -q 2i + 2 — i n contrast to my result © — , even 
after having derived the correct integral ©. N. Kaiser claimed |96j . that 
the 1/a term in the ERE could be neglected, from wich they obtained in [2E] 
an Q cm -dependence of the cross section close to threshold. An alternative 
explanation for the mysterious threshold energy dependence of modern cal- 
culations hidden in computer codes might be another systematic mistake, 
i.e. the use of a 2-body phasespace instead of a 3-body phasespace. E.g. 
J. A. Niskanen states in [21], that a ~ (a 2 p/)/(l + (p/ a) 2 ), where "pf in the 
numerator is the momentum dependence of the phase space". Well, from 
the arguments given above we know, that the momentum dependence must 
be p 2 ! Finally I want to mention here a semantic problem in the definition of 
%: in [SI] H.O. Meyer et al. state correctly, that r] n is "the largest possible 
center-of-mass pion momentum (with nucleons at rest relative to each other) 



divided by the pion mass" , which yields rj^ := y X(s, m 2 , s™ n )/(4 s m 2 ). In 
Horowitz et al. (1994) [02] it is stated, that q' is not k, but "the pion momen- 
tum in the final pp center-of-mass system" , which would yield the definition 



^Horowitz ._ ^/A(s,m 2 ,s 1 mi ")/(4sf in m 2 ) = r] w yjs/{2m p ) 2 . Even if both 
quantities vanish at threshold, they increasingly differ at higher energies. 
7^ = 0.557 yields e.g. r]^ orowitz = 0.603. 



38 Even a infinitesimal value for 1/a would yield o~(pp — > ppn°) oc i]% close to threshold! 

39 This is quite surprising, as J. Niskanen states in |93| . that H.O. Meyer (see e.g. p. 
657 in 1851). G.A. Miller, P.U. Sauer |7fci| are aware, that the cross section "dependence at 
threshold is not a to t = arf . . . " . In their case this seems to be due to the energy depen- 
dence of their specific form of (see e.g. 14911 ) and not due to a corrected phasespace. 
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The Miller- Sauer problem 

The energy dependence of the total cross section for pp — ► ppn derived 
by G.A. Miller and P.U. Sauer (1991) [75] — the standard calculation for 
direct meson production — is not reproducable by means discussed in their 
publication. The fall-off of the total cross section for rj n > 0.4 is much too 
steep compared to Effective Range or Coulomb Modified Effective Range 
calculations or for a result obtained by a calculation with wavefunctions 
as solutions of Schrodinger equations with strong and Coulomb potential. 
The direct meson production is only possible through FSI and ISI and is 
therefore a crucial test for a correct description of FSI and ISI. 

The problem of wrong application of the two potential formalism 

Many people nowadays describe the influence of Coulomb interactions in FSI 
by the CFS-CMERE [SH [2H1 listed in © and derived from a two potential 
formalism. I.e. the total S-wave phaseshift 5o(k) is given by 5q(k) = <to(k) + 
5q(k), while co(re) is the S-wave Coulomb phaseshift. In the construction of 
Tpsi we have — according to the two potential formalism — to use 5q(k) 
and not <5 c (re). Therefore Tpsi ls given in the OWMA according to (|34|) by 

h + hr c ^ 2 -^-2K V H(r ] )-i K C H V ) + ) ' 

(46) 

The limit cxo(re) — > performed by the majority of experimental and theo- 
retical researchers yields results, which one would get, if one does not apply 
the two potential formalism properly. As discussed above the presence of 
Coulomb phaseshifts o~o(k) introduces essential singularities in the phase- 
space integrations, when calculating cross sections, which can only be re- 
solved, if the Coulomb potential is regularised consistently. 

The numerator or Moalem et al. problem 

Going back to the ideas of K. Brueckner et al. j^j, i.e. that Tpsi °^ VK r = 0) 
or equivalently Tpsi — {Nexp(i 5) sin 6)/k = N/(k (cot 6 — i)), it has been 
very early clear, that the denominator of enhancement factors is determined 
by the ERE, i.e. by k (cot 5 — i) = —a -1 + ... — i re. The unknown numerator 
of Tpsi has been fixed either by hand, i.e. by setting the dimensional con- 
stant N to some suitable value, or — to be free of ambiguities — Tpsi was 



Tfsi - 
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chosen to be the inverse Jost function 40 , i.e. Tpsi = 1//o(k). Even if the 
inverse Jost function yields — apart from the correct denominator of Tpsi 
- a K-dependent numerator, it can be observed in literature, that many 
authors approximate the numerator of Tpsi (after performing a Taylor ex- 
pansion in k) by the term, which is O(k ), i.e. by a constant, while keeping 
the full K-dependence of the denominator (see e.g. [TSl 1571 1581 155] r 1 . 

Among others A. Moalem et al. (1996) @U HUU] and R. Shyam (1999) 
|42| use within an Effective Range Approach complex scattering length and 
effective range parameters, Moalem for the np system in the final state, 
Shyam for the K + A system in the final state. Both authors use a Coulomb 
modified CFS formula |321 129j to describe Coulomb interactions and both 
authors observe mysterious "oscillations" of the cross section as a function 
of energy close to threshold. How can this be, especially as both use different 
methods for enhancement factors: Moalem uses Tpsi = JV//c(cot 8q (k) — i), 
while Shyam uses the inverse Jost function Tpsi = l/fo(n). The problem 
is, that both authors don't use the correct numerator 42 in the enhancement 
factor, which is determined by the OWMA, i.e. they must use Tpsi = 
Kcot 5q{k) / {n{cot 5o(k) — i)) = (Re/o(^))//o(^)) which stabilizes the en- 
ergy dependence very close to threshold. Both authors don't apply the two 
potential formalism correctly by skipping the Coulomb-phaseshifts and the 
Rutherford amplitude. On the other hand one can observe, that — due 
to the different numerators of Tpsi — the penetration factor C 2 (rj) enters 
differently in the "more correct" (above) and wrong (below) expression: 



KCot£g(«) - ± + \ r c k 2 + ... -2KnH(rj) 



(47) 



«(cot5§(/0 - i) - i + I r c k 2 + ... -2KrjH{rj) - inC^{-q) 

N ^ NC 2 (r,) 

K(cot<5§(«) -i) ~^ - i + \r c k 2 + ... -2KrjH(rj) - inC^n) 

The Titov et al. problem 

As a side effect of the purpose of their calculation Titov et al. (2000) |101j re- 
considered to some extend the Inverse Jost Function Method for the descrip- 
tion of ISI and FSI. The problem arising in their work is, that they use for 

40 This choice can be found in nearly every textbook of scattering theory (see e.g. 012])- 
It should be noted, that this nontrivial approximation changes the energy dependence 
of the total cross section at larger energies. 

42 By their choice both authors set the numerator of the enhancement factor to a con- 
stant: Moalem replaces KCot<5o(K) by a constant, Shyam performs first the replacement 
Re/o(fi:) — > 1, then he removes the k 2 term of the ERE in the numerator ending up with 
an enhancement factor like Moalem, i.e. Tpsi = iV'//t(cot <5o( K ) — *)• 
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the description of the pp FSI the effective range parameters a pp = —7.8098 
fm and r pp = 2.767 fm, even not using a CMERE. In order to describe their 
inverse Jost function they essentially use a simple nuclear "shape indepen- 
dent" ERE Kcot5(K) = —a -1 + r k? /2. In such a nuclear ERE without 
Coulomb they should use for the pp system the well known effective range 
parameters a pp = —17.1 fm and r pp ~ 2.8 fm (see e.g. p. 254 in |1U2| . p. 425 
in |43| . p. 115 in 30 ). The relation beween nuclear and Coulomb modified 
effective range parameters is discussed e.g. in |103| or on p. 39 in |104j . 



The wavenumber problem 

The threshold production of heavy mesons in NN collisions involves high rel- 
ative energies in the initial state and high energy- and momentum-transfers 
at the meson-production vertex. In quantitative calculations it is therefore 
important to define quantities carefully and such, that they are valid also at 
high energies within a relativistic framework. This is, why it is important to 
mention, that standard EREs are expansions in the Lab-wavenumber and 
not in momenta 43 . 



The phaseshift problem 

Phaseshifts are described in the Lab-system! This is not only true for pro- 
ton nucleus scattering, but also for pp scattering, where the cm-system 
seems more convenient. If phaseshifts are theoretically determined from 
a Schrodinger equation (e.g. for the pp-system (p. 419ff in [IS]), it makes a 
large difference at higher energies, whether one uses Lab- or cm-quantities: 



1 d ( r 2 d \ ! K 2 L{L + 1) arn p \ r> , u 
r 2 dr \ dr 1 r 2 r 



= rn p J2VLS,L'S'(r)R^ Sf ( K ^)- ( 48 ) 

US' 

if cm = r, fLab = r /2). Le. use Lab-quantities K = k and T-^^ = uj p (2k) — 
m p — 2 k 2 /m p (Final State (FS)) or K = kq and = lu p (2kq) — m p ~ 

43 It is therefore not correct, if on p. R241 in |78| it is stated, that "fc is the relative baryon 
momentum". The authors should have written, that k is the relative Lab-wavenumber, 
which makes of course only at high relative energies a difference! On page R242 in |78| 
the authors define the Gamov factor by C(q) = 2 n 7 3 /(exp(2 n j q ) — 1) with f q — /-ta/q 
and state, that "fi and q are the reduced mass and the relative momentum of the two- 
body subsystem in which the Coulomb interaction occurs . . ." . This definition is also not 
exact at higher energies! They should have defined 7 9 = a/v with v being the relative 
Lab-velocity of the two-body subsystem in which the Coulomb interaction occurs. 
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2 Kq /m p (Initial State (IS)) and not the cm-quantities K = 2m p K/ 'y/s\ and 
T = yfs\ — 2m p (FS) or K = 2m p kq/^s and T = — 2m p (IS)! 



The transition operator problem 

In traditional threshold production models for direct 7r°-production in pp — > 
pp7T° close to threshold people don't use one consistent 7r°-production oper- 
ator. There are e.g.: 

oo 

o 

oo 

Tfi oc Jdrr 2 (R^(K,r)rj (K n or)^ + ^R^(K ,r)^m 



o 

oo 



Tfi oc «vo J drr 2 (R^(K,r))* n^or) R$( K(h r) -^M ■ (49) 



which are more or less related by a partial integration 44 . For a correct 
description of FSI we need the correct transition operator. 



The data reduction problem 

For a FSI reduction of their data experimentalists tend to multiply enhance- 
ment factors describing the FSI of possible particle subsystems in the final 
state. I quote here e.g. J. Smyrski et al. |1U5| . who states: "The ppr\ FSI 
can be factorized into pp{f pp ) and prj(f pn ) factors and integrated over the 
available phase-space volume p 3 : a(e) ~ / f P p(q P p) • f PV (q PlV ) ■ fprifaptv) d /°3, 
. . ." . For the moment there is no obvious theoretical reason to choose a prod- 
uct of final state factors rather than a sum or anything else to describe the 
FSI! Well, the leading terms of a Faddeev expansion would suggest a the 
coherent square of a sum. It is an important theoretical issue in the future, 
to give a satisfactory answer to this still open problem. 



44 It is unclear, why various terms appearing during this partial integration should vanish 
or be neglected. The transition operator describes mainly the FSI in the NN-system, as 
join^o r) representing the produced 7r°-meson is nearly constant close to threshold. 
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8 Conclusions and outlook 



Theoretical aspects 

We come to the conclusion, that — in order to be able to relate quanti- 
tatively and conclusively theoretical quantities to experimental data — one 
either has has to calculate the process NN — ► NNX close to threshold com- 
pletely, e.g. within a three-body Faddeev approach, or there is a need of a 
quantitative theory to separate and estimate Initial-/Final-State effects. It 
is unsatisfactory and demanding, that — in spite of the vast related activ- 
ities of theoreticians during the passed about 60 years — such a theory for 
estimating Initial- /Final- State effects does not exist yet! This theory also 
has to be extended to differential cross sections for processes involving unpo- 
larized and/or polarized particles close to threshold to extract independent 
information from data. Even being under control in the differential descrip- 
tion of elastic scattering of charged particles the application of the theory 
of Coulomb scattering in initial and final states of particle number changing 
reactions is — due to the infinite range of the Coulomb interaction and the 
offshellness of the particles in the interaction zone — a highly nontrivial 
problem, requiring well defined regularization schemes, especially close to 
the particle production threshold, at which produced charged particles have 
a long time to interact. The development of an adequate theoretical — if 
possible, relativistic — formalism is an outstanding task. In order to handle 
interference effects and inelasticities a consistent field theoretical treatment 
of intermediate resonant states is demanded. An adequate formalism for 
the treatment of Fermionic |15| llf)8[ I109| and Bosonic |110| resonances has 
been recently proposed. Keeping in mind the high relative energies and 
momenta of the incoming nucleons and the high momentum transfers in- 
volved in production processes of heavy mesons, relativistic formulations 
for nucleon induced meson production close to threshold are crucial. This 
includes the consideration of boost effects and e.g. tensor components, neg- 
ative energy and parity components of the deuteron wavefunction within 
quantitative calculations. For the consistent relativistic treatment of bound 
of quasi-bound composite systems in the initial and final state the author 
refers to the new available formalism called ABSSM 73 j. It is a pitty to 
observe, that the majority of theoretical models stays nonrelativistic, espe- 
cially because most of the traditional models have been developed to describe 
pion production only. Let's go for solving relativistic Faddeev problems! For 
the iViV 7r-system I refer e.g. to |1U61 11U7] . Not everything what is written 
should be believed! 
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Experimental aspects 



Thanks to experimental facilities like the COSY in Jiilich, the CELSIUS 
ring in Uppsala or the IUCF in Bloomington the accuracy of cross section 
data related to nucleon induced meson production processes close to thresh- 
old improved significantly within the last ten years. It is a challenging, 
yet unexpected difficult task for theory to handle the experimental results 
quantitatively. In order to obtain as much independent experimental in- 
formation as possible to constrain theory it is highly desirable to continue 
with the experimental investigation of differential cross sections, polarization 
and polarization transfer observables, production of heavy mesonic systems 
with and without strangeness, processes involving a deuteron beam or tar- 
get and/or further light and heavy nuclei. For a quantitative estimate of 
the effect of ISI on nucleon induced meson production reactions close to 
threshold an experimental determination of NN scattering phaseshifts and 
inelasticities at high relative energies is crucial. As proposed in it is 
highly desirable to compare experimentally uncharged and charged reaction 
channels like e.g. nn — > nnX and pp — ► ppX to learn something about 
the influence of Coulomb effects in the Initial/Final State and the effective 
range parameters of uncharged strong interacting particle systems. In the 
meantime we ask for being patient with slowly theoreticians! 45 
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